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ABSTRACT 

We present sensitive and high angular resolution CO(l-O) data obtained by the Combined Array 
for Research in Millimeter-wave Astronomy (CARMA) observations toward the nearby grand-design 
spiral galaxy M 51. The angular resolution of 0.7" corresponds to 30 pc, which is similar to the typical 
size of Giant Molecular Clouds (GMCs), and the sensitivity is also high enough to detect typical 
GMCs. Within the 1' field of view centered on a spiral arm, a number of GMC-scale structures are 
detected as clumps. However, only a few clumps are found to be associated with each Giant Molecular 
Association (GMA), and more than 90% of the total flux is resolved out in our data. Considering the 
high sensitivity and resolution of our data, these results indicate that GMAs are not mere confusion of 
GMCs but plausibly smooth structures. In addition, we have found that the most massive clumps are 
located downstream of the spiral arm, which suggests that they are at a later stage of molecular cloud 
evolution across the arm and plausibly are cores of GMAs. By comparing with Ha and Paa images, 
most of these cores are found to have nearby star forming regions. We thus propose an evolutionary 
scenario for the interstellar medium, in which smaller molecular clouds collide to form smooth GMAs 
at spiral arm regions and then star formation is triggered in the GMA cores. Our new CO data have 
revealed the internal structure of GMAs at GMC scales, finding the most massive substructures on 
the downstream side of the arm in close association with the brightest H II regions. 
Subject headings: galaxies: individual (M 51 or NGC 5194) - galaxies: spiral - ISM: clouds - ISM: 
molecules 



1. INTRODUCTION 

As all the stars are formed from the interstellar 
medium (ISM), the evolution of physical and chemical 
conditions of ISM is an important clue to understand 
star formation, its history, and galaxy formation. Spi- 
ral galaxies have been studied often with this goal; they 
are actively forming stars in spiral arms and the se- 
quence of the ISM evolution from molecular gas to young 
stars can be investigated across the arms. There have 
been many discussions on the nature of the spiral struc- 
ture. iLin fc Shul (|1964l ) described it as a quasi-stationary 
structure sustained by density waves, while tidal waves 
due to interactions with a companion galaxy or shear mo- 
tions in the gal actic disk can also create transient spiral 
structures fe.g.. iToomre k, Toomrdfl97^ . 

With its proximity, nearly face-on disk, and grand- 
design spiral structure, the Whirlpool galaxy, M 51 or 
NGC 5194, has been a popular target for both obser- 
vational and theoretical studies. Str ong streaming mo- 
tions found in th e two spiral arms (|Aalto et all Il999t 
iShettv et al.l2007T> indicate the existence of density waves 
in this galaxy, while it has been suggested that the spiral 
structure is not fully stable across tim e or radius (e.g., 
iMeidt et all 120081 : iDobbs et all 120101 ) . An interaction 
with the companion galaxy, NGC 5195, has been thought 
to enhance or even trigger the prominent spiral structure 
in the disk of M 51. Molecular gas in this grand-design 
spiral galaxy traced by dust lanes and CO emission lines 
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has been fou nd upstream of star forming regions trace d 
by Ha (e.g., IVogel et all 119881 : IRand fc Kul Wnil fl990T) . 
iTilanus et alF i 19881 ) found similar spatial displacement 
between non-thermal and thermal radio continuum; spi- 
ral arms seen in the non-thermal component are coin- 
cident with dust lanes while most of the thermal com- 
ponent peaks agree with H II regions. These offsets, 
also found in other spiral galaxies, are consistent with 
the density wave theory and represent s the timescal e 
for molecular clouds to form stars (e.g., lRobertsl fl969). 
Variations of physical properties across the arms thus 
corr espond to the evol ution of the ISM in spiral galax- 
ies. lEgusa et all ((2009) utilized offsets between CO and 
Ha to derive a star formation timescale and pattern 
speed in nearby spiral galaxies. 

Given that the typical width of spiral arms is about 
500 pc, investigating the ISM evolution across the arms 
requires images with a high spatial resolution, which 
has been difficult especially for radio observations. In- 
stead, comparison between arm and interarm regions 
has been used to diagnose how the spiral structure af- 
fects the ISM properties. From the 12 CO(1-0) data 
covering the entire disk of M 51 taken with the Com- 
bined Array for Research in Millimeter- wave Astronomy 
CAR MA) and Nobeyama 45m telescope, iKoda et all 
2009, hereafter K09) found that Giant Molecular Asso- 
ciations (GMAs; > 10 7 M ; IVogel et all 1 19881) are only 
in arm regions while Giant Molecular Clouds (GMCs; 
10 5 ~ 6 M ; ISolomon et al1ll979ah are in both arm and 
interarm regions. On the other hand, bright young 
stars delineate narrow spiral arms while older stars are 
distributed m ore smoot hly over a galactic disk (e.g., 
ISandagdll96l . For M 51. IScoville eFall (pOOl cataloged 
H II regions from a high-resolution Ha image taken with 
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the Hubble Space Telescope (HST) and found that about 
half of them are confined to the arm region which occu- 
pies only quarter of the disk area. It is thus expected 
that GMAs play a critical role in massive star formation 
in the spiral arms. 

However, the typical size of GMAs and H II regions 
are quite different; > 200 pc and ~ 30 pc, respectively. 
Thus, the key to understanding the initiation of massive 
OB star clusters is the internal structure of GMAs. K09 
measured the surface density of GMAs to be similar to 
that of typical GMCs and suggested that GMAs are not 
simply overlapping of GMCs but are distinct structures 
where the entire volume is virtually filled with molecu- 
lar gas. Nevertheless, the spatial resolution of their data 
(~ 200 pc) is not high enough to resolve GMAs. Since 
previous CO observations with high spatial resolution 
similar to the size of GMCs and H II regions have only 
been done for the central region of galaxies or for dwarf 
and flocculent spiral galaxies in the Local Group, the in- 
ternal structure of GMAs in grand-design spiral galaxies 
has not been investigated before. 

To resolve GMAs and detect GMC-scale structures 
(~ 40 pc) in spiral arms, we have carried out new 
CARMA observations in the CO lines toward M 51 with 
a very high angular resolution (~ 0.7"). In this paper, 
we present results from our new CO observations and 
discuss the internal structure of GMAs and the relation- 
ship of the brightest substructures to the spiral structure 
and star for ming regions. At the a dopted distance of M 
51 (8.4 Mpc: lFeldmeier et al.lfl997b . 1" corresponds to 41 
pc. 

2. OBSERVATIONS AND DATA REDUCTION 

The CO (1-0) line observations were carried out at 
CARMA in November through December 2008. This 
array has six 10 m dishes and nine 6 m dishes and thus 
provides 105 baselines. The full width at half maximum 
(FWHM) of the primary beam is about 60" and 100" 
for the 10 m and 6 m dishes, respectively. The B-array 
configuration, whose typical angular resolution is 0.8" at 
3 mm, was selected to resolve GMAs and detect GMC- 
scale structures in M 51. A single pointing was centered 
on the brightest GMA in the K09 data located about 30" 
from the galactic center in the south-west direction. The 
pointing center coordinate is (RA, DEC)=(13:29:50.16, 
+47:11:28.39) in J2000. We used the hybrid correlator 
mode (three 62 MHz bands for the source and three 500 
MHz bands for the calibrator) in order to achieve high 
spectral resolution for the source and high sensitivity for 
the calibrator. Channel width for the source is 0.977 
MHz or 2.54 km s _1 at the frequency of 115 GHz. The in- 
termediate frequency (IF) was set to include the 12 CO(l- 
0) line in the upper side band (USB) and 13 CO(1-0) in 
the lower side band (LSB) . Data were taken through the 
standard cal-source cycle with 5 minutes integration on 
the phase calibrator, 1153+495, which is 16 degrees away 
from the source, and 15 minutes on source. The pass- 
band observations were done at the beginning of each 
track with a brighter source, and the data were used 
to calibrate phase and amplitude offsets between narrow 
and wide bands as well as the normal passband calibra- 
tion. The total observing time is about 34 hours with 
five tracks. In terms of the weather condition, the phase 
rms was ~ 150-200 /im and the opacity at 230 GHz was 



~ 0.3, which is typical for 3 mm observations with the 
B-array configuration. 

Data reduction was done with the Multichannel Im- 
age Recons truction, Image A nalysis and Display software 
(MIRIAD; iSault et al.lll995h . The flux of the phase cal- 
ibrator was derived to be 0.81 Jy at 114.75 GHz (USB) 
and 0.74 Jy at 110.02 GHz (LSB) from the compari- 
son with bright quasars whose flux had been derived 
from other observations with planets. The uncertainty 
in the flux calibration is estimated to be 15% (~ 0.1 
Jy). The synthesized beam from the natural weighting 
is 0.73" x 0.60" for USB and 0.77" x 0.62" for LSB, corre- 
sponding to ~ 31 pcx25 pc at the distance of M 51. The 
noise rms (a) measured at central 40" x 40" in emission- 
free channels of the dirty map is 12 mJy beam -1 or 2.5 
K for USB and 7.2 mJy beam" 1 or 1.5 K for LSB, with 
the channel width of 5.08 km s _1 as Hanning smoothing 
was applied. The baseline length ranges 23.9 kA to 362 
kA, which corresponds to 8.6" to 0.60" in angular resolu- 
tion or 360 pc to 24 pc. Given that the typical width of 
molecular spiral arm seen in the K09 CO data is ~ 10" 
or 400 pc, most of the structures smaller than the arm 
down to 24 pc should be detected. 

The dirty map was CLEANed by the use of the 
MIRIAD task mossdi with the flux cutoff of la. In or- 
der to shorten the time needed to CLEAN, we applied a 
mask made by the combination of (i) the field of view, 
and (ii) mask from the K09 cube data. First, we created 
the sensitivity map from the primary beam patterns by 
the task mossen and defined the field of view as where the 
sensitivity is lower than 1.5 times the one at the pointing 
center. The diameter of this field of view is 53.25" and 
pixels outside have been masked out to avoid confusion 
with elevated noise. We made the RA-DEC, RA-VEL, 
and DEC-VEL projections of the K09 cube and manu- 
ally drew a polygon around emissions to make the mask 
for the entire galaxy. Within the field of view, only the 
region inside the K09 mask has been used for CLEAN- 
ing. We refer to this CLEANed dataset as B-array data 
hereafter. 

The integrated intensity or moment [Iqo — 
J Tcodv) map was created by integrating the B-array 
data cube with the velocity range of 385.5-563.2 km s _1 . 
The flux cutoff for making the Iqo map was set to be 
proportional to the local sensitivity; la at the pointing 
center and 1.5<7 at the edge of the field of view. The 
resulting Iqo map for the 12 CO(1-0) data is shown in 
the left panel of Figure [IJ while the B-array field of view 
is indicated as a box on the K09 Zoo map in the right 
panel. 

3. RESULTS AND DISCUSSION 

3.1. Detecting GMC-scale structures 

In order to identify each mol ecular cloud structure , the 
clumpfind method provided by Willi ams et al.l (|1994f ) has 
been applied to the CLEANed 12 CO (1-0) data cube with 
the lowest level and contouring step of la. Both of these 
two parameters are constant within the field of view. We 
have then defined detections as clumps (i) within the field 
of view, (ii) with the peak flux whose local signal-to-noise 
ratio (S/N) > 4, and (iii) with the FWHM in velocity 
larger than 0.5 pixel (or 2.54 km s -1 ). We have calcu- 
lated this local S/N, considering the variable sensitivity 
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Figure 1. Left: Integrated intensity map of 12 CO(1-0) data from the B-array CARMA observations. The unit of the gray scale is Jy 
beam -1 km s — 1 . The synthesized beam (0.73" X 0.60") is shown as the ellipse on the bottom left corner. Circles superimposed on the 
gray-scale image indicate clumps detected by the clumpfind, wit h the line width and color corresponding to the local S/N of the peak flux; 
thin blue is for S/N > 4 and thick red is for S/N > 5. See £|3.1l for detail. Right: The field of view of the left image is shown as the box 
overlaid on the lower resolution 12 CO(l— 0) integrated intensity map from K09. Note that the color scale for the two images is not the 
same. 



within the field of view. 

The selected 34 clumps are indicated by circles in the 
left panel of Figure Q] superimposed on the moment 
map. The line width and color represents the local S/N of 
each clump; thin blue is for 5 > S/N > 4 while thick red 
is for S/N > 5. The radius measured by the clumpfind is 
0.5"-1.4" or 10-60 pc, corresponding to the radius of the 
circles in the figure. The FWHM in RA and DEC also fall 
in the same range with the radius and about one-third 
of the selected clumps are resolved. While most of the 
clumps are found in the arm, larger and more massive 
clumps are located only on the west side of the spiral 
arm. Assuming a trailing arm, this corresponds to the 
downstream side. On the other hand, some of smaller 
clumps are found in the upstream interarm region (east 
side of the field of view). In H'S.21 we discuss it in more 
detail by comparison with the K09 data. 

From Figure [TJ we have noticed that there are several 
emission peaks not selected as clumps. Most of them are 
close to the edge of the field of view and thus are pre- 
sumably elevated noise, while some appear sequentially 
on upstream side of the spiral arm. We have examined 
all the clumps with S/N > 3 and found that they are 
superposition of several weak- and narrow-line clumps 
at different velocities. As some of them seem to be as- 
sociated with H II regions (see Figure [5]), they perhaps 
are real molecular clouds, but we need a higher sensitiv- 
ity to confirm that and do not include them for further 
discussion. 

For the 13 CO(1-0) line in the LSB data, the same 
procedures as the 12 CO(1-0) data have been applied 
and three clumps with S/N > 4 have been detected. 



None of them, however, are found to be associated with 
the 12 CO(1-0) clumps. Lowering the S/N cutoff to 3, 
we have detected 162 clumps and found two clumps 
with S/N ~ 3.6 located near 12 CO(1-0) clumps; one 
is located close to the pointing center and coincides 
with the largest and most massive 12 CO(1-0) clump 
(marked as 'd' in Figure^, while the other is located 
close to the two 12 CO(1-0) clumps at ~ 10" north of 
the pointing center ('b' and 'c' in Figure A ra- 
tio of the peak temperatures of corresponding clumps 
is T peak ( 13 CO)/T pcak ( 12 CO) ~ 0.4, while a ratio of the 
total flux is J F 1 tot ( 13 CO)/i ;l tot ( 12 CO) - 0.2. The aver- 
age integrated i ntensit y ratio for the Galaxy is 1/5.5 
([Sob mon et all Il979bf ) . implying that these massive 
clumps detected in M 51 are in similar condition to 
GMCs in the Galaxy. As these 13 CO (1-0) structures do 
not meet the S/N criterion of for the clumps, we discuss 
only the 12 CO(1-0) results hereafter. 

The conversion factor from the 12 CO(1-0) integrated 
intensity to the H2 column density (iV(H2) = -^co-^co) 
has been measured in a number of ways. For the Galaxy, 
Xqo — 3 x 10 20 cm -2 (K km s -1 ) -1 was derived by com- 
parin g the virial mass and the 12 CO(1-0) flux of GMCs 
(e.g- JSolomon et al.lll987l;lYoung fc Scovilldl 19911 ). More 
recently, smaller values (~ 1.8 x 10 20 cm~ 2 (K km 



s 1 ) 1 ) have been re ported by comparison with the 7-ray 



or far infrared data ( 


Strong & MattoxHl996t IDame et al.l 


120011 iGrenier et all 


20051 IHeveretal.1 120091) derived 



GMC masses from the 13 CO (1-0) emission lines under 
the local thermodynamic equilibrium (LTE) assumption 
and found that the LTE masses are factor of ~ 5 smaller 
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than the virial masses. The inconsistencies in Xqo values 
are due to the combination of uncertainties and/or inva- 
lidities of the models and assumptions introduced and 
imply difficulties in the Xqo determination. 

As for the Galaxy, derived Xqo values for M 51 depend 
largely on data an d methods. While t he vir ial masses of 
GMAs derived by IRand fc Kulkarffil {1990) agree with 
X gp = 3 x 1 20 cm " 2 (K km s -1 ) -1 , those derived 
by lAdler et"aT1 (|1992l ) based on a different CO dataset 
suggest X C o = 1-2 x 10 20 cm" 2 ( K km s" 1 )" 1 . Even 
smaller values have been suggested dGarcia-Burillo et al.1 
1991 IR.andl 119931: INakai fc Kund llQQfih. More recently. 
Schinnerer et al.l ()2010[ ) applied the large velocity gradi- 
ent (LVG) modeling to 2" resolution multi-transition CO 
data and derived X C o = (1-3 - 1.9) x 10 20 cm" 2 (K km 
s -1 ) -1 . In this paper, we adopt Xqo = 1-5 x 10 20 cm~ 2 
(K km s^ 1 )" 1 . The variations of Xqo within the field 
of view as a function of galactocentric radius should be 
small enough to be neglected, as the data only cover 53" 
or 2.2 kpc. 

The noise rms (la) in the 12 CO(1-0) data thus corre- 
sponds to A/h 2 = 2.5 x 10 4 Mq within the synthesized 
beam (~ 30 pc) and one velocity channel (5.08 km s _1 ). 
Though the absolute mass of molecular gas could be 
changed by a factor of two or more according to adopted 
Xqo values, the sensitivity in mass is still smaller than 
typical GMCs in any case. Furthermore, the conclusion 
of this paper does not depend on the Xqo value, since 
we only focus on the spatial and kincmatical information 
and discuss relative values of mass from the B-array and 
K09 data. 

3.2. Relationship to larger structures 

In the left panel of Figure [3 contours from the K09 Ico 
map are overlaid on the B-array Iqo image. We have ap- 
plied the clumpfind algorithm to the K09 data cube and 
defined GMAs as clumps with M > 10 7 M . Following 
the selection criteria for the B-array data, clumps on the 
data edge have been excluded. These selected GMAs are 
shown as open black circles in the right panel of Figure [2] 
together with the detected clumps in the B-array data in 
the same manner as Figure [TJ The background image is 
the K09 moment map. From this Figure, we find that 
GMAs are not filled with the GMC-scale clumps. As 
the la sensitivity of the B-array data is about 10 times 
smaller than the typical mass of GMCs, we should have 
detected tens or more of them if GMAs were confusion 
of GMCs. This result strongly indicates that GMAs are 
not confusion nor overlapping of GMCs. 

The total flux within the field of view of the B-array 
data has been measured to be 67 Jy km s _1 and 1.1 x 10 3 
Jy km s _1 for the B-array and K09 data, respectively. As 
the K09 data include the single-dish observations with 
Nobeyama 45m telescope, it should represent the total 
12 CO(1-0) flux. The missing flux in the B-array data 
is thus 94%, suggesting that most of the CO flux in the 
spiral arm of M 51 should come from either extended 
structures or molecular clouds smaller than GMCs and 
supporting the aforementioned indication that GMAs are 
not a collection of GMCs. 

Another important clue from comparison of the two 
Ico maps is that their peak positions do not coincide; 
massive chimps detected in the B-array data are shifted a 



few arcseconds toward the downstream side of the global 
spiral arm seen in the K09 map. To investigate this dif- 
ference further with the velocity information, we have 
made spectra of B-array and K09 data at the positions 
of the selected clumps. The box size for the average spec- 
tra was set to be the same as the diameter of each clump 
derived by the clumpfind. In Figure^ sample spectra for 
massive clumps are displayed. The flux difference men- 
tioned above is clearly seen also in the spectra. From 
these spectra, we have found that massive clumps have 
offsets in peak velocity between the two datasets. We 
have defined the velocity offset as 

V ois = Ki(B) - V peak {K09), (1) 

where V^B) is the peak velocity for clumps measured by 
the clumpfind and V r pca k(K09) is the peak velocity in the 
K09 spectra. The V f s is plotted against the mass of each 
clump in Figure [4] It is clear from this plot that almost 
all of the clumps with Mh 2 > 5 x 10 5 Mq have negative 
V f s . Considering the galactic rotation, negative veloci- 
ties correspond to the downstream side. Therefore, mas- 
sive clumps are kinematically as well as spatially located 
downstream of the spiral arm. This difference indicates 
that the massive clumps are in a later evolutionary stage 
compared to the larger scale structures or GMAs. We 
refer to these massive clumps as GMA cores hereafter. 

Smaller clumps are on the other hand located preferen- 
tially on upstream of the spiral arm, and little emission is 
found in the interarm region on downstream. This distri- 
bution of small clumps indicates that they are somehow 
destroyed through the evolution across the spiral arm. 
Before forming stars, the small clumps (M ~ 10 5 Mq) 
could collide with each other to form larger clouds such 
as GMAs. Meanwhile, they would be dissociated into 
atomic gas or broken up to even smaller clouds (M < 10 5 
Mq) after forming stars. 

3.3. Relationship to star forming regions 

As already mentioned in SJTJ star forming regions are 
located on the downstream side of the spiral arms. Since 
we have found that GMA cores are also on downstream, 
it is worth comparing these two components. 

The Ha data taken under the HST program (ID: 10452, 
PI: S. Beckwith) were obtained via the HST website 4 
and cover ~ 7' x 10' inclu ding the entire galact ic disk 
and the companion galaxy (M utchler et al.l 120 05). The 
Paa image taken for another HS T program (ID : 7237, 
PI: N. Scoville) was presented by iScoville et ail (|2001[ ) 
and covers central 3' of the disk. The astrometry of the 
Paa image was found to be slightly off compared to the 
Ha image. In order to match their coordinates, we mea- 
sured the peak position of bright H II regions seen in both 
Ha and Paa images, and adjusted the coordinates of the 
Paa image by the use of the task ccmap of the Image 
Reduction and Analysis Facility (IRAF). 5 Although the 
fitting rms from ccmap is as small as 0.17" and 0.13" for 
RA and DEC, respectively, the errors in the coordinates 
could be larger as H II regions are extended especially in 

4 http:/ /archive. stsci.edu/prepds/m51/ 

5 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 2. Left: Contours from the K09 moment map overlaid on the B-array moment map. The synthesized beams for these two 
datasets are indicated by ellipses on the bottom left corner. Right: Clumps from the B-array data (as circle s wi th b lue or red line) and 
GMAs from the K09 data (as circles with black lines) superimposed on the K09 moment map. See i|3.1l and 13.21 for their definition. 
Clumps whose spectra are shown in Figure[3]are marked with 'a'— 'd'. 
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Figure 3. Example of 12 CO(1-0) spectra for selected massive 
clumps, which are indicated in the right panel of Figure [2] Thin 
red lines indicate the K09 data while black thick lines indicate the 
B-array data. Note the difference in flux scale between the two 
datasets. 

the Ha image and extinction gradients might affect the 
flux distribution and thus the peak positions. We esti- 
mate the effective uncertainty in the coordinates to be 
about 40 pc or 1", which is the typical size of GMCs and 
H II regions. In Figure [5j circles indicating the selected 
clumps are superimposed on the Ha and Paa images. 

As expected, bright H II regions are located on the 
downstream side of the CO spiral arm seen in the lower 
resolution K09 data. GMA cores found in the high res- 
olution B-array data, on the other hand, coincide with 
them. This correlation indicates that GMA cores are the 
site of massive star formation, which is consistent with 
the result in the previous subsection that they are in a 
later stage of the molecular gas evolution. Meanwhile, 
we have found that not all the H II regions are accom- 
panied by the cores, which implies that GMA cores are 
short-lived or are not prerequisites for massive star for- 



Figure 4. Plot of V D f s defined in Eq. JTJ v.s. mass of the selected 
clumps. Larger points represent clumps with S/N > 5. A few 
datapoints outside the V f s range are not shown. 

mation. 

Given that only the 12 CO(1-0) data are available 
at this resolution, physical conditions such as density 
and temperature of GMA cores cannot be determined 
uniquely; high density and high temperature can both 
produce high I CO - The largest core close to the point- 
ing center, (marked as 'd' in Figure [5]), gives a clue on 
this issue. From Figure [SJ we have found that it does 
not have a counterpart in Ha but in Paa, indicating a 
high extinction and thus a high column density. From 
the mass (2.5 x 10 6 M@) and radius (55 pc) measured by 
the clumpfind, the average H 2 column density is derived 
to be 2V(H2) = 1.6 x 10 22 cm~ 2 for this core. Assuming 
A V /N K = 5.3 x 1(T 22 mag cm 2 ijBohlin et al.lll97cl . the 
estimated extinction is as high as Ay = 17 mag. The 
most plausible detection of 13 CO(1-0) at this location is 
also consistent with the high density. Therefore, at least 
for this core, we conclude that the high Iqo represents 
high column density of molecular gas. This agreement 
in high extinction and density also confirms that the co- 
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Figure 5. Selected clumps presented as circles (same as Figure[T} on the HST He* (left) and Paa (right) image. One-third from the right 
edge of the field of view, where no clumps have been detected, is not displayed. 



ordinate adjustment for the Paa image is correct within 
the core size or ~ 1". 

Additionally, it is clear from Figure [S] that the distri- 
bution of H II regions is quite similar in the Ha and 
Paa images, indicating only a few regions (such as the 
most massive core mentioned above) are significantly af- 
fected by extinction. This means that Ha is a good tracer 
of the locations of current star formation in most cases 
and is preferable over extinction-free 24/im for measur- 
ing offsets between molecular arm and young stellar arm 
(|Egusa et al.ll2009T ) as it generally provides higher angu- 
lar resolution. 

3.4. Evolution scenario across a spiral arm 

The results presented above reveal the structure and 
distribution of molecular clouds at 30 pc scale for the 
first time in grand-design spiral galaxies. Based on 
their properties and relationship to the global (kpc) spi- 
ral structure, we here propose a scenario for the ISM 
evolution across the spiral arm. We should note here 
that smaller spatial scale across the arm corresponds to 
shorter timescale, which enables us to describe the sce- 
nario in unprecedented detail. 

On the upstream side of the spiral arm, molecular 
clouds of Mh 2 ~ 10 5 M Q collide with each other and 
probably with smaller clouds to grow into GMAs (Mh 2 > 
10 7 Mq). GMAs are presumably distinct and smooth 
structures and a few GMC-scale cores (~ 10 6 Mq) are 
formed within each GMA, which turn to be the site of 
OB star formation as they evolve. If not all the gas in 
cores would turn into stars, remaining gas in the core 
should break up to smaller molecular clouds (< 10 5 Mq) 
or become dissociated into atomic hydrogen by strong 
radiation from y oung stars. From the star formation 
law presented by iKennicutt et all ([2007) for M 51, the 



gas consumption timescale measured as the ratio of the 
total gas mass to the star formation rate is ~ 10 9 yr. 
Sin ce the star formatio n timescale is much shorter (~ 10 7 
yr; lEgusa et al.l I2009D than this timescale, most of the 
gas in GMAs is not used by the star formation activity. 
Such gas remains molecular through the evolutionary se- 
quence, considering that the molecular fraction to the 
total hydrogen has been estimated as 70-80% even in 
the interarm regions (K09). 

4. SUMMARY 

We have carried out high resolution CO observations 
toward a spiral arm in the nearby spiral galaxy M 51 
using the radio interferometer CARMA, in order to (i) 
resolve GMAs, (ii) detect GMC-scale structures, and (iii) 
understand their relationship to the global spiral struc- 
ture. For the 12 CO(1-0) line, the angular resolution and 
noise rms (la) are 0.7" x 0.6" (or ~ 30 pc) and 12 mJy 
beam -1 (or 2.5 K) with the velocity width of 5.1 km s _1 , 
respectively. Assuming the CO-to-H2 conversion factor, 
4cr corresponds to 10 5 Mq within the beam, which is 
smaller than the mass of typical GMCs. The spatial res- 
olution is similar to the typical size of GMCs, and thus 
the data quality is high enough to detect CMC-like struc- 
tures. 

Within the 1' field of view, a number of GMC-scale 
structures referred as clumps are detected in the 12 CO(l- 
0) data. For the 13 CO(1-0) line, only two marginal de- 
tections have been identified. By comparing with lower 
resolution 12 CO(1~0) data (K09), we have located only 
a few clumps in each GMA and derived missing flux to 
be larger than 90%. These results indicate that GMAs 
are not mere confusion of GMCs but are plausibly single 
and smooth structures. Among the detected clumps, we 
have found that the most massive (> 5 x 10 5 Mq) clumps 
are located at the downstream side of the spiral arm spa- 
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tially as well as kinematically. This displacement indi- 
cates that they are at a later stage of the ISM evolution 
across the spiral arm and plausibly are cores of GMAs. 
Spatial coincidence between the massive clumps and H II 
regions supports this hypothesis. Smaller clumps on the 
other hand are found preferentially in the upstream side 
of the spiral arm. 

We thus propose an evolutionary scenario for the ISM 
in spiral galaxies; as smaller clouds in interarm re- 
gions approach to spiral arms, they cluster or collide to 
form GMAs, which are smooth and discrete structures. 
Within a GMA, a few massive and GMC-scale substruc- 
tures referred to as cores are formed at a later stage of 
its lifetime. Star formation is triggered within the GMA 
cores, and finally the remaining gas in the core is dissoci- 
ated into atomic gas or broken up to even smaller clouds 
after stars are formed. 

The sensitive and high resolution CO data from our 
new CARMA observations have revealed the internal 
structure of GMAs at GMC scale and its relationship 
to the spiral arm structure for the first time in exter- 
nal galaxies. This result is an important step to fully 
understand the ISM evolution scenario in spiral galaxies. 
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us a draft of her paper prior to publication. Support for 
CARMA construction was derived from the states of Cal- 
ifornia, Illinois, and Maryland, the James S. McDonnell 
Foundation, the Gordon and Betty Moore Foundation, 
the Kenneth T. and Eileen L. Norris Foundation, the 
University of Chicago, the Associates of the California 
Institute of Technology, and the National Science Foun- 
dation. Ongoing CARMA development and operations 
are supported by the National Science Foundation un- 
der a cooperative agreement, and by the CARMA part- 
ner universities. This research is partially supported by 
HST-AR-11261.01. 
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